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We consider quantum particle or Gaussian polymer confinement between two surfaces and in 
cylinders with sinusoidal undulations. In terms of the variational method, we show that the quan- 
tum mechanical wave equations have lower ground state energy in these geometries under long 
wavelength undulations, where bulges are formed and waves are localized in the bulges. It turns out 
correspondingly that Gaussian polymer chains in undulated boxes or tubes acquire higher entropy 
than in exactly flat or straight ones. These phenomena are explained by the uncertainty principle 
for quantum particles, and by a polymer confinement rule for Gaussian polymers. If membrane 
boxes or tubes are flexible, polymer-induced undulation instability is suggested. We find that the 
wavelength of undulations at the threshold of instability for a membrane box is almost twice the 
distance between two walls of the box. Surprisingly we find that the instability for tubes begins 
with a shorter wavelength compared to the "Rayleigh" area-minimizing instability. 



I. INTRODUCTION 

Recently much attention has been paid to the struc- 
ture and dynamics of polymer chains restricted by two 
surfaces, or restricted in cylindrical pores These 
conditions are relevant to a broad class of applications 
and biological functions, such as filtration, gel perme- 
ation chromatography, heterogeneous catalysis, oil recu- 
peration. 

When a polymer is confined to a smaller space, the 
entropy loss of a chain is higher. Here we call this simple 
rule the polymer confinement rule. It has been employed 
casually in the literature, however, it has a close connec- 
tion with the uncertainty principle in quantum mechan- 
ics. Indeed, on a theoretical side, quantum mechanical 
wave equations have been applied to elucidate the behav- 
ior of the Gaussian chains when the radius of gyration of 
a chain is much larger than these structures. Therefore, 
the solutions of wave equations in some geometries are 
useful both for mesoscopic quantum physics and poly- 
mer physics. 

The most prominent recent example is the calcula- 
tion for curved geometries. Goldstone and Jaffe || have 
shown that the bend of two and three dimensional tubes 
with a constant cross section lowers the ground state en- 
ergy of a quantum particle constrained in the tubes. Cor- 
respondingly, the cylindrical bend of two parallel walls 
(2D tube) with a constant width reduces the entropy loss 
Yaman et al. have shown in their series of works 
that entropic interactions between curved membranes 
and macromolecules such as flexible chain molecules and 
rigid rods may change the bare elastic constants of single 
membranes and bilayers 0. 

In this paper, we focus on tubes in two and three di- 



mensions with sinusoidal undulations keeping tube vol- 
umes constant. This lowers the ground state energy of 
quantum particles, and thus raises the entropy of Gaus- 
sian polymer chains. Corrugated walls induce additional 
kinetic energy along the walls, however, if bulges are 
formed, the waves are less confined in the transverse di- 
rection and localized in the bulges, which may decrease 
the total energy because of the uncertainty principle. In 
the same way, according to the polymer confinement rule, 
polymer chains may favor the undulations of tubes. 

In section II, wc first describe the polymer confinement 
rule, an intuitive argument for why a polymer in a con- 
fined space loses entropy. This approach is microscopic, 
which is different from the well-known scaling argument 
Second, we provide the calculation procedure for 
Gaussian chains using quantum mechanics which is given 
m Ref. @. In section III, we elucidate the effect of undu- 
lations using the variational method. We show that the 
wave equations have lower energy under long wavelength 
undulations. Correspondingly, long Gaussian chains ac- 
quire higher entropy with undulating constraints than 
with exactly flat or straight ones. For the three dimen- 
sional case, we compare the undulation effect with the 
Rayleigh area- minimizing instability In the final 

section, we discuss implications of polymer-mediated en- 
tropic force for membrane boxes and tubes. 



II. CORRESPONDENCE TO QUANTUM 
MECHANICS 



1 



A. The polymer confinement rule 



then, the problem reduces to the wave equation, 



The entropy reduction of a confined random walk of N 
steps between two parallel flat walls (distance Aa; ) can 
be derived intuitively through the following microscopic 
reasoning. The number of steps that span the distance 

1/2 

between the walls is described by Aa; ~ IN , and the 
gyration radius in a free space R g ( 3> Aa;) is expressed 
as R g ~ IN 1 / 2 , where I is a step length. Here ~ implies 
that the numerical factors have been ignored. 

Steps touching the walls should return. This require- 
ment is the source of the entropy reduction |)j . The chain 
reflects off the walls about N/N' times. Let z be the num- 
ber of nearest positions on a lattice for a step. Then the 
number W of total configurations of the confined chain 
is represented by 



W ~ z N - N l N '(z/2) N l N ' = z N (l/2) N / N \ 



(1) 



where the combination factor has been omitted. Then 
the entropy loss due to the walls for the Gaussian chain 
is AS = keN/N' hi2~ k B R 2 J{Ax) 2 . We write 



(Aa;) 2 A5 ~ k B R 2 



(2) 



which corresponds to the uncertainty principle: 

(Aa;) 2 A.E ~ (Aa;) 2 (Ap) 2 /m ~ h 2 /m, (3) 

where E, x, p, and m are the energy, position, momen- 
tum and mass of a quantum particle, and h is the Planck 
constant. On this basis, one can say that a polymer 
tends to escape from narrower spaces into wider or open 
spaces; in other words, the polymer tends to localize in 
wider spaces of the confined geometries, the same as with 
quantum particles. 

B. Analytical theory for Gaussian chains 

For a random walk, we consider the total number of 
paths that connect r and r' with N steps, z N G(r' , r, N), 
where z is the number of neighboring sites. The bound- 
ary condition is 



G(r',r,0) - V, r 



(4) 



It is easy to show G(r', r, N) is the solution of the diffu- 
sion equation: 



^-^)^'.r..V, = ,.. 



where I is the lattice constant, or the step length. 
With the eigenfunction expansion: 



(5) 



G(r', r, N) = I 3 ^ *;(r')*„(r) exp 



Nl 2 Ej, 



(6) 



V 2 * n (r) + E n V n (r) = 



(7) 



with the boundary conditions. As is known, the eigen- 
functions satisfy the orthogonality and completeness con- 
ditions: 



(8) 
(9) 



dr#^(r)# n (r) = S mn , 



Since the continuous limit is 



Um5 r i r /l 3 = S(r-r'), 



Eq.(|6|) satisfies the boundary condition Eq.(|y). 
For excited states (i > 1), when an equality 



(10) 



Rl (Ei — Ex) » 1, 



(11) 



is satisfied, the situation is called ground state domi- 
nance, which is our interest in this paper. Then we have, 



G(r,r',JV) w*5(r')*o(r)exp - 



Nl 2 Ex 
~6 



(12) 



The partition function is a sum over all configurations: 
Z = y J drdr'G{r,r',N). (13) 

Integrating out all uninteresting degrees of freedom, 
and we obtain the main term in the associated entropy 
change, 



S = /cr In Z —k] 



Nl 2 E x 



(14) 



For a Gaussian chain between two flat walls, the dis- 
crete part of the eigenvalues of Eq. (Q) is 



E n = T-r , 2 , 71= 1,2,3,- •• 

(Aa;) z 



(15) 



To obtain the main term, we ignore the continuous part 
of the eigenvalues associated with two directions along 
walls. Then Eq.dTl]) is fulfilled. Thus, we have 



S w —h 



7T 2 Rl 



6 (Ax) 2 ' 
This is consistent with Eq.(||). 



(16) 
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III. VARIATIONAL PROOFS 

Let C be a tube in two dimensions (box) or in three 
dimensions. We will consider the wave equation, 



(V 2 + E)i>(r) = 0, 



(17) 



in C subjected to the Dirichlet condition on walls: 
ij}(dC) — 0. The following calculations are applicable 
to both quantum particles and Gaussian polymers. 
We define 



(18) 



where D is the dimension of the system. Our aim is to 
show 



AM = a[>p] - <J > 0, 



(19) 



where ctq corresponds to the ground state for the exactly 
flat or straight case. Since perturbations are sinusoidal, 
the integrals can be limited to a single wavelength, that 
of the perturbations. We compute a critical wavelength, 
down to which A[ip] is positive. Our evaluation gives 
an upper bound for the critical wavelength, because we 
employ a trial function for -0 in Eq.(|l8|). 



A. In two dimensions 

As shown in Fig.l, a box is defined between two par- 
allel lines described by 



We will prove, with the variational method, that long- 
wave undulation lowers the total energy. Choose a can- 
didate function with a variational constant 5: 

ib(x,y) =cos(—) (1 + Sesw.ax), 0<5<l/e. (22) 
v w i 

S > implies localization at bulges. The denominator of 
Eq.(pl)is 



<l.r I dyij? = -y 



where 

7 = 1 + (s + 
It is elementary to calculate 

2 



e 2 > 0. 



dy 2 
d 2 ^ 
dx 2 



w ■ f 7T V\ ^ w ■ C K V\ 
7T— ysml — -2?r— ysml — 

w \ w J w a \ w J 



y cos I — 

V w 



(1 + Se sin ax) 



zoeaTT — J?/ sin ( — ) cos ax 



w 

ea" cos 



5'" 2 



\ w ) 



\ w I 
sin ax. 



Then, a brief calculation yields 

1 f ^ ^p 2 ^ 71-2 (1 + ^esinax) 2 
d J dy 2 2d 2 1 + e sin ax 



(23) 



(24) 



iw) C0S (~^) ( 1 + Ses[nax )' ( 25 ) 



(26) 



(27) 



y±(x) = ±-(l + esina:r), 0<e<l (20) 



with a width undulation function w(x): 
w(x) = d(l + e sin ax), 



(21) 



where a = 2tt/X, X is the period of the undulation, and 
d is the mean distance. Note that the volume is kept 
constant. 



d/2 


y 

o" x 










-6/2 







FIG. 1. Undulated box in two dimensions: A is the period 
of the undulation, and d is the mean distance between walls. 



1 f , .d 2 ip 7r 2 + 6 e 2 a 2 cos 2 ax . . , 2 

d J = " 24 l + e S inax {1 + S€SmaX) 



{ + S j e 2 a 2 sin 2 ax 
S 2 e 3 a 2 



(28) 



2 

3S 2 e 3 a 2 



cos ax + 



sin 3 ax 



cos ax sin ax 

2 

(ea 2 + 2d"ea 2 ) 

- sin ax. 



4 4 
Using formulae given in Appendix A, we obtain 

7*„M = - J (-71=1 - 2SA + 62 A ) > 

ja x [ip] = —a 2 BD - 25a 2 ( ( — - CD 



$V -+CD 



(29) 
(30) 



where a[tp] = o y [ip] + ^[-0], and where 



3 



.4 



B = 1- yl 

/2 



C = B 



D = 



2 

•6 



e 2 

- 1 = — 4 
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For e > 0, we evaluate the deviation A[ip] 
-w 2 /d 2 . To show A[^] is positive, we write as 



(31) 
(32) 
(33) 
(34) 
for ctq = 



jA[ip] = 



26 



7T 
d 2 

ze! 
rf 2 



^(-40D + e 2 

7T 2 

— (40D + 2e 2 ) 



(35) 



From this, it easy to show that there exists 8 that makes 
A[ip] positive in the long- wave undulation limit: A — > oo; 
8 = 1 for instance. Therefore, the ground state energy 
should be lower than that in the flat plane case. No- 
tice that as is seen in Eq.(29) the decrease in energy 28 A 
appears in the transverse direction. 

In order to estimate the critical wavelength Ao, an ap- 
proximation up to order e 2 is 



( 7 /e 2 )AM 



-2<5 2 
'2d 2 



TV 

A 2 

- 2D 



28 



7T~ 
tf 2 



7T 

A 2 



A 2 ' 



(36) 



Maximizing A[tp] by changing the variational constant 
5, we compute the critical wavelength from A[tp] = 0. 
Thus, 



Ao 
d 



4tt 2 
3 
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2.014. 



(37) 



The critical wavelength is about two times the width of 
the box, implying that the bulge is at least the size of the 
width between walls. It is a quite reasonable value, be- 
cause, in terms of quantum mechanics, if the undulation 
is shorter than Ao, then the kinetic energy along walls 
will be higher than the confinement energy. 



B. In three dimensions 

It was shown by Plateau, later pursued by Rayleigh, 
then known as the Rayleigh instability, that the undula- 
tion of a tube with wavelength exceeding its circumfer- 
ences reduces the surface area. Therefore, it is interest- 
ing to compare the polymer-mediated instability and the 
Rayleigh instability. 



Consider a tube of unperturbed radius Rq and an un- 
dulation function (Fig. 2): 



R(z) — R c (l + esinaz), < e < 1, 



(38) 



where a = 2tt/\, and R c is determined by the constant 
volume condition of the tube. Thus. 



R 2 c = Rl{l + e 2 /2)-\ 



(39) 




FIG. 2. Undulated tube in three dimensions: A is the pe- 
riod of the undulation. 

In order to remind readers of the Rayleigh instability, 
we calculate the surface area of the tube given by 



A(e; A) =2rrJ Rdl 



(40) 



We expand the square root for small perturbations and 
use Eq.(p9|). The area is 



A(e; A) w 27ri? A 



(2tt.Ro) 2 
A 2 



1 



(41) 



Therefore, if the wavelength is greater than the circum- 
ference of the tube, the area is less than that of the un- 
perturbed perfect cylinder. 

Now, we will prove that a long-wave undulation lowers 
the total energy with the variational method. Choose a 
trial function with a variational constant 8: 

^( r , z ) = J (l + Sesinaz), 0<8<l/e. (42) 

where Jo(r) is the 0-th Bessel function: 

^+^Ml) 2 MfH. <«) 

and where a = 2.40483 is the smallest zero of Jq(x). The 
denominator of Eq. JTs| ) is 



dz I dr2-Krip~ z = tt\R z c (3 7, 
'0 



where j3 = J 2 (a) and 



1 = l + \{8 2 +AS+l)e 2 + l8 2 e A . 



(44) 



(45) 
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It is elementary to calculate 



IV. DISCUSSION 



d 2 ip 1 dip 
dr 2 r dr 

d 2 ^ 
dz 2 



(S) 2j °(l^) (1 + ,5eSillaz) ' (46) 



R" i (ar\ R' 2 . / ar 



(1 + <5esinaz) 



cos az 



(47) 



Sea 2 , Jo sinaz 



where, J'q{x) — dJo(x)/dx, R' = dR(z)/dz, and R" = 
d 2 R{z)/dz 2 . 

Using the calculation in Appendix B, we then obtain 



2 C 2\ 



a 2 (2 + S 2 e 

w 

4 + a 2 



(48) 



o . 2 6 2 e 2 13 + a 2 r2 4 

e 2 + <5e 2 H 1 5V 

6 2 24 

For e > 0, we evaluate A[ip] for cr = —o?jR\ = -b 2 

7 AM 



-6 2 e 2 o/l 13 + a 2 



8 



+ (2b 2 - a 2 )5 - a 



2 24 

2 4 + « 2 
6 



An approximation up to e 2 is 



^<5 2 + (26 2 -a 2 )(5-a 24 + a2 



G 



(49) 



(50) 



Maximizing A[ip] by varying 6, we calculate the critical 
wavelength from A[tp] = 0. Thus, we have 



where 



Aa 4 x 4 - Aa 2 x 2 - = 0, 



A 



2ttRo aa ' 



(51) 



(52) 



is the ratio of the critical wavelength to that of the 
Rayleigh instability. Then the minimum xq > is 



•i'O 



4+a 2 



2a 2 



0.493. 



(53) 



The critical wavelength is shorter than that of the 
Rayleigh instability. 



We have shown that quantum mechanical particles 
confined in undulated boxes or tubes have lower energy 
when the wavelength is greater than certain values, com- 
parable to the width between walls or the radius of tubes. 
We have explained that the effect can easily be inter- 
preted by the uncertainty principle. Quantum mechan- 
ical calculations immediately imply that long Gaussian 
chains in undulated boxes or tubes acquire higher en- 
tropy than in exactly flat or straight ones. Furthermore, 
it can be explained by the polymer confinement rule es- 
tablished in this paper, which is quite analogous to the 
uncertainty principle. 

This polymer-mediated entropic force may play an im- 
portant role in a number of systems: polymers in cell 
membranes, vesicles, microemulsions, and polymers con- 
fined in lamellar or cylindrical phases of surfactant and 
homopolymers blending into those of block copolymer 
systems. For instance, deformable flat membrane boxes 
with constant width d containing the Gaussian chains 
(R g 3> d) are unstable against the undulations. In the 
same way, cylindrical tubes with radius R containing 
the Gaussian chains (R g 3> R) are unstable against the 
undulations. From Eqs. (p6|) and (|50|), it is easy to see 
that the entropy gain increases with increasing undula- 
tion amplitude. Therefore, the undulation amplitude will 
grow until it is balanced by the elastic restoring force 
of the membrane, or until the membrane ruptures. In 
addition, the critical wavelength is shorter than that of 
the Rayleigh area-minimizing instability; the polymer- 
mediated interaction may trigger the undulation insta- 
bility. 

When the bare elastic modulus of a membrane is not 
small, for long wavelength undulations, the effect softens 
the surface tension of the membrane box. On the other 
hand, for short wavelength undulations, the effect hard- 
ens the surface tension because of the reduction of the 
entropy of confined polymers. It is quite remarkable that 
the discrimination point of softening or hardening (in- 
stabilizing or stabilizing) is determined by the distance 
between membranes. 

One should be careful to interpret our results. First, 
because of the volume preserving condition, the undu- 
lation is not a simple expansion of membrane boxes or 
tubes by the thermal motion of polymers. Second, the 
calculation of wave equations does not correspond to the 
case when the wave length of undulations is longer than 
confined polymer sizes in boxes or tubes. 

The polymer-mediated entropic force without undu- 
lations is proportional to d -3 , where d is the distance 
between walls, which decays slower than Van der Waals 
attractive interactions. This force is something like the 
excluded volume interactions of multimembranes systems 
known as the Helflich repulsive interaction. Finally, the 



5 



polymer confinement rule for excluded volume chains is 
modified as (Ax) 5 / 3 AS ~ k B R 5 g /3 . We expect that the 
same entropic effect will exist for excluded volume chains. 



APPENDIX A 

To tackle the integration of Eqs.(|27|)-(|28|), we define an 
integral: 



< m >= 



2tt J_2l 1 + esinax 
By using a formula, 

dx 2 



1 + esinx ^/l — e 2 



■. tan 



_ x / tan | + e 

Wi - 



we have 



(54) 



(55) 



drr 3 J 2 {cr) = i[r 4 J 2 (cr) + (r 4 - l r 2 )J 2 ( C r) 

+ -r 3 J (cr)J 1 (cr)]. (64) 
c 



Then, we obtain 



'drrVof^Ur' 



2 ' 



\RJ U \RJ 

"*^(=)-^(i-£>- 



n- 



(65) 
(66) 
(67) 



< f >= 



Hence, it yields 



< sin arc > = —A/e, 

< sin 2 ai > = A/e 2 , 

< cos 2 ax > = B/e 2 , 
< cos 2 ax sin ax > = — C/e 3 , 

< cos 2 ax sin 2 ax > = C/e 4 , 



APPENDIX B 



where A, i? and C are defined in Eqs.(BfJ) — (33) 



(56) 



(57) 
(58) 
(59) 
(60) 
(61) 



In this appendix, we provide some integrals of the 
Bessel functions. By using Jq(x) = — Jx(x) and 
(xJi(x))' = xJq(x) 1 it is easy to verify the following in- 
tegrals: 



aYrJoV) = [- (J 2 (cr) + J 2 (cr))], (62) 
J drr 2 Mcr)4(cr) = h^iV)], (63) 
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